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Abstract. Disproportionation in 3d compounds can give rise to an unconventional 
electron-hole Bose liquid with a very rich phase diagram from a Bose metal, charge 
ordering insulator to an inhomogeneous Bose superfluid. Optimal conditions for the 
disproportionation driven high-T c superconductivity are shown to realize only for 
several Jahn- Teller d n configurations that permit the formation of well defined local 
composite bosons. These are the high-spin d 4 , low-spin d 7 , and d 9 configurations 
given the octahedral crystal field, and the d , high-spin d e configurations given the 
tetrahedral crystal field. The disproportionation reaction has a peculiar "anti-Jahn- 
Teller" character lifting the bare orbital degeneracy. Superconductivity in the d 4 and d 6 
systems at variance with d 1 , d 7 , and d 9 systems implies an unavoidable coexistence of 
the spin-triplet composite bosons and a magnetic lattice. We argue that unconventional 
high-T c superconductivity observed in quasi-2D cuprates with tetragonally distorted 
CuOg octahedra and iron-based layered pnictides/chalcogenides with tetrahedrally 
coordinated Fe 2+ ions can be a key argument supporting the disproportionation 
scenario is at work in these compounds. 
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1. Introduction 

The origin of high-T c superconductivity [T] is presently still a matter of great controversy. 
Copper oxides start out life as insulators in contrast with BCS superconductors 
being conventional metals. Unconventional behavior of these materials under charge 
doping, in particular, a remarkable interplay of charge, lattice, orbital, and spin 
degrees of freedom, strongly differs from that of ordinary metals and merely resembles 
that of a doped semiconductor. Novel non-copper based layered high-T c materials 
such as iron oxyarsenide LaOFeAs [2] reveal normal and superconducting state 
properties very different from that of standard electron-phonon coupled " conventional" 
superconductors. 

We believe that unconventional behavior of high-T c superconductors can be 
consistently explained in frames of a so called dielectric scenario EJ [5] that implies 
an instability of the respective parent compound with regard to d-d charge transfer 
(CT) fluctuations with formation of local composite bosons. This implies realisation 
of the Ogg-Schafroth's real-space pairing model for the Bose-Einstein superconducting 
condensation [6j [7j. At first sight the d-d disproportionation 

3d n + 3d n ->• 3d n+1 + 3d" 1 - 1 (1) 

whose energy is usually estimated to be U (Udd), mean intra-atomic electron-electron 
(Mott-Hubbard) repulsion energy, is rather uncommon and a priori much more expensive 
in narrow-band 3d transition-metal compounds as compared with compounds containing 
broad-band 6s 1 ions such as Bi 4+ , Pb 3+ or Tl 2+ due to their strong tendency to form 
6s° and 6s 2 closed shells [8]. Experimental values of the minimal energies of the d-d CT, 
or Mott-Hubbard transitions in different 3d oxides as usually derived from the optical 
gap measurements, A^f ~ 2-4 eV (see, e.g., Ref. [9] and references therein) yield effective 
IPs which are sizeably less than typical U^d ~ 5-10 eV, however, far from the so-called 
"negative-t/" regime [8]. Yet A^ is a minimal energy cost of the optically excited d-d 
disproportionation or electron-hole formation due to a direct Franck-Condon (FC) CT 
transition. The question arises, what is the energy cost for the thermal excitation of 
such a local d-d disproportionation? The answer implies first of all the knowledge of 
the relaxation energy, or the energy gain due to the lattice polarization by the localized 
charges that can exceed leV [10]. For instance, in insulating cuprates such as La2Cu04 
the true (thermal) d-d charge transfer gap appears to be as small as 0.4-0.5 eV rather 
than 1.5-2.0 eV as derived from the optical gap measurements [5J. In other words, parent 
cuprates should be addressed to be d-d CT unstable systems. 

At present, the CT instability with regard to the d-d disproportionation and 
formation of "negative-f/" centers [8] is believed to be a rather typical property for a 
number of perovskite 3d transition- metal oxides such as CaFe0 3 , SrFe0 3 , RNi0 3 [TT||T2]. 
furthermore, in solid state chemistry one consider tens of disproportionated systems [13J. 

Classical disproportionation seems to be observed in the perovskite ferrate CaFe0 3 , 
where the average formal valence of the iron ion is Fe 4+ (3d 4 ) is known to show a gradual 
charge disproportionation 2Fe 4+ — )-Fe 3+ +Fe 5+ below 290 K with antiferromagnetism 
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below 115 K according to Mossbauer studies [13] • The single magnetic hyperfine pattern 
for isostructural ferrate SrFeC>3 at 4K, on the other hand, indicates seemingly a 
rapid electron exchange between Fe 3+ and Fe 5+ ions, and the hyperfine field coincides 
approximately with the average value of the corresponding parameters for CaFe0 3 . 
Probably, in this 3d 4 oxide we meet with a manifestation of quantum effects and 
formation of a quantum dimerized phase such as a "valence bond solid" composed of S- 
or P-type neutral dimers described by wave functions: 



respectively. The Fe Mossbauer measurements point to a charge disproportionation 
in Sr 3 Fe207, which has a double-layered perovskite structure [T5] . 

Rare earth nickelates RMO3 (R is a trivalent rare earth ion) exhibit a first order 
metal-insulator phase transition (MIT) upon cooling with distinct signatures of the 
charge disproportionation [121 [16] ■ At T^r < Tmit these exhibit a further transition to a 
long-range antiferromagnetic ordered state with a nonzero magnetic moment for one of 
the two nonequivalent Ni sites [17] . A clear fingerprint of charge disproportionation is 
the breathing-type distortion of metal-oxygen octahedra since the different charge states 
of the transition-metal ion take different ionic radii (see, e.g., Ref. [T2]). 

However, a large body of 3d materials does not manifest the charge instability 
so clearly as (Ca,Sr)Fe03 or RN1O3. For instance, the perovskite manganites RMnC>3 
with the same as in ferrates 3c? 4 configuration of Mn 3+ ions reveal a "hidden" CT 
instability [T8l [T9] . At first glance the disproportionation in manganese compounds 
is hardly possible since manganese atom does not manifest a "valence-skipping" 
phenomenon as, e.g. bismuth atom which can be found as Bi 3+ or Bi 5+ , but not 
Bi 4+ , with a generic bismuth oxide BaBiC>3 to be a well-known example of a charge 
disproportionated system (see, e.g., Ref. [20] and references therein). However, strictly 
speaking, sometimes manganese reveals a valence preference, e.g., while both Mn 2+ and 
Mn 4+ are observed in MgO:Mn and CaO:Mn, the Mn 3+ center is missing [2Tj. 

The reason for valence skipping or valence preference observed for many elements 
still remains a mistery. Recently, Harrison [22] has argued that most likely traditional 
lattice relaxation effects, rather than any intra-atomic mechanisms [8] (specific behavior 
of ionization energies, stability of closed shells, strong screening of the high-charged 
states) are a driving force for disproportionation with formation of "negative-f/" centers. 

Earlier it was argued that the disproportionated system can form an unconventional 
electron-hole Bose liquid whose phase diagram incorporates different phase states from 
classical, or chemical disproportionated state to quantum states, in particular, to 
the unconventional Bose-superfluid (superconducting) state [3] ■ Regrettably physicists 
have paid remarkably little attention to the problem of valence disproportionation and 
negative-f/ approaches ("chemical" route!) which are surely being grossly neglected in 
all present formal theoretical treatments of HTSC Speaking about a close relation 
between disproportionation and superconductivity it is worth noting the text-book 
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example of BaBiC>3 system where we unexpectedly deal with the disproportionated 
Ba 3+ + Ba 5+ ground state instead of the conventional lattice of Ba 4+ cations [8]. The 
bismuthate can be converted to a superconductor by a nonisovalent substitution such 
as in Ba 2+ i_ x K + a; Bi03. At present, this system seems to be the only one where 
the unconventional superconductivity is related anyhow with the disproportionation 
reaction. 

In the paper we concern ourselves with great points of a so-called " disproportiona- 
tion" scenario in 3d compounds which was addressed earlier by many authors, however, 
by now it was not properly developed as its consistent description in terms of conven- 
tional Hubbard type models is complicated, if any. Our final goal is to ascertain the 
criteria for a motivated search of the 3d systems which are the most promising ones 
as high-temperature superconductors. The paper is organized as follows. In Sec. II ad- 
dressing a simple model of the CT unstable 2D system we demonstrate a potential of 
the disproportionation as a driving force for superconductivity. In Sec. Ill we analyze 
the features of the disproportionation reaction for different 3d n pairs and do find 3d n 
configurations optimal for the disproportionation driven superconductivity. In Sec. IV 
we address the electron structure of familiar 3d compounds with the optimal 3d n config- 
urations and point to high-T c superconductivity observed in quasi-2D cuprates and iron 
pnictides/chalcogenides to be a result of the d-d disproportionation. A short conclusion 
is made in Sec.V. 

2. Simple toy model of the mixed valence system and the 
disproportionation driven superconductivity 

2.1. Pseudospin description of the model mixed-valent system 

Valent electronic states in strongly correlated 3d oxides manifest both significant 
correlations and p-d covalency with a distinct trend to localisation of many-electron 
configurations formed by antibonding Me 3d-0 2p hybridized molecular orbitals. The 
localisation effects are particularly clear featured in the crystal field d-d transitions 
whose spectra just weakly vary from diluted to concentrated 3d oxides. An optimal 
way to describe valent electronic states in strongly correlated 3d oxides is provided by 
quantum-chemical techniques such as the ligand field theory |23j which implies a crystal 
composed of a system of small 3d-cation-anion clusters. Naturally, such an approach 
has a number of shortcomings, nevertheless, this provides a clear physical picture of 
the complex electronic structure and the energy spectrum, as well as the possibility 
of a quantitative modelling. In a certain sense the cluster calculations might provide a 
better description of the overall electronic structure of insulating 3d oxides than different 
band structure calculations, mainly due to a better account for correlation effects and 
electron-lattice coupling. 

In the 3d oxides with the disproportionation instability (JTJ we should consider 
as a minimum three different many-electron configurations, or different valent states, 
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with different types of the charge transfer. One strategy to deal with such a mixed- 
valence system is to create model pseudospin Hamiltonians which can reasonably well 
reproduce both the ground state and important low-energy excitations of the full 
problem. Standard pseudospin formalism represents a variant of the equivalent operators 
technique widely known in different physical problems from classical and quantum lattice 
gases, binary alloys, (anti)ferroelectrics,.. to neural networks. The formalism starts with 
a finite basis set for a lattice site (triplet of M°,M ± centers in our model, see below). 
Such an approach differs from well-known pseudospin-particle transformations akin 
Jordan- Wigner [24J or Holstein-Primakoff [25] transformation which establish a strict 
linkage between pseudospin operators and the creation/annihilation operators of the 
Fermi or Bose type. The pseudospin formalism for electron systems generally proceeds 
with a truncated basis and does not imply a strict relation to fermion operators that 
obey the fermionic anti-commutation rules. 

To demonstrate the perspectives of the disproportionation in driving the 
superconductivity we address hereafter a simplified toy model of a mixed-valence 
system with three possible stable valence states of a cation-anion cluster (Cu0 4 , Mn0 6 , 
FeAs4,...), hereafter M: M°,M ± , forming the charge (isospin) triplet and neglect all 
other degrees of freedom saving only the quantum charge one. Similarly to the neutral- 
to- ionic electronic-structural transformation in organic charge-transfer crystals (see, e.g., 
paper by T. Luty in Ref. [26]) the system of charge triplets can be described in frames 
of the S=l pseudo-spin formalism [1] . To this end we associate three different valence 
charge states of the M-center: M°, M with three components of the S = 1 pseudo-spin 
(isospin) triplet with M s = 0, +1,— 1, respectively. Having in mind quasi-2D cuprates 
we associate M° ) M ± centers with three charge states of the CUO4 plaquette: a bare 
center M°=Cu04~, a hole center M + =CuO|~, and an electron center M~=Cu04~, 
respectively. 

However, the physically simple toy mixed-valence system should be described by 
rather complex effective pseudospin Hamiltonian as follows [I] 

H = £(A,4 - h { S iz ) + v {j s iz s jz + 

i i<j 

Y^iDffiSi+Sj- + 5i_5 i+ ) + D%\T i+ Ti„ + T^T j+ )} 

i<j 

+ J2Usls]_ + sls]^ (2) 

i<j 

with a charge density constraint: J2i(Si Z ) = An, where An is the deviation from 
a half-filling (N M + = N M ~). Two first single-site terms describe the effects of a bare 
pseudo-spin splitting, or the local energy of M 0,± centers. The second term may be 
related to a pseudo- magnetic field hj || Z, in particular, a real electric field which acts 
as a chemical potential. The third term describes the effects of the short- and long- 
range inter-site density- density interactions including screened Coulomb and covalent 
couplings. The last three kinetic energy terms in ([2]) describe the one- and two-particle 



Perspectives of the disproportionation driven superconductivity 



6 



hopping, respectively (T± = {S z , S±}). All the parameters have a clear physical meaning. 
The energy Act for the creation of uncoupled electron and hole centers, or effective 
parameter U<id is: Act = 2 A, while the energy Aeh for the creation of coupled electron 
and hole centers, or EH-dimer, is: A EH = 2 A — V nn . 

One should note that despite many simplifications, the effective pseudospin 
Hamiltonian fl2]) is rather complex, and represents one of the most general forms of the 
anisotropic S — 1 non-Heisenberg Hamiltonians. Its real spin counterpart corresponds 
to an anisotropic S=l magnet with a single ion (on-site) and two-ion (bilinear and 
biquadratic) anisotropy in an external magnetic field. Spin Hamiltonian (J5J) describes 
an interplay of the Zeeman, single-ion and two-ion anisotropic terms giving rise to a 
competition of an (ant i) ferromagnetic order along Z-axis with an in-plane XY magnetic 
order. Simplified versions of anisotropic S = 1 non-Heisenberg Hamiltonians have been 
investigated rather extensively in recent years pTl [28] . 

Our system is characterised by several order parameters. These are two classical 
(diagonal) order parameters: (S z ) being a valence, or charge density with an electro- 
neutrality constraint, and (Sf) being the density of polar centers M ± , or "ionicity". 
In addition, there are two unconventional off-diagonal order parameters (S+) and (T + ) 
related to two different types of a correlated single-particle transport, and one order 
parameter related to a two-particle transport. It is worth noting that the 

operator creates an on-site hole pair, or composite boson, with a kinematic constraint 
= 0, that underlines its "hard-core" nature. For real S=l spin systems we deal with 
conventional "linear" magnetic order parameters (S z ) and («S+) related to diagonal (Z-) 
and off-diagonal (XY-) orderings and several spin-quadrupole order parameters: (S 2 ), 
(T + ), and (S+). Recently, the spin-quadrupole (spin-nematic) order (S^) has been 
investigated in Ref. [28J. 




The last three terms in (j2J) representing the one- and two-particle hopping, 
respectively, are of primary importance for the transport properties, and deserve special 
attention. Two types of one-particle hopping are governed by two transfer integrals 
D( 1,2 \ respectively. The transfer integral tL = {Dif + D^) specifies the probability 
amplitude for a local disproportionation, or the eh-pair creation: 

M° + M° — >■ M ± + M T , 

and the inverse process of the eh-pair recombination: 

M ± + M T -> M° + M°, 

while the transfer integral t"j = (Dij — D^) specifies the probability amplitude for a 
polar center transfer: 



or the motion of the electron (hole) center in the lattice of M° -centers or motion of the 
M°-center in the lattice of M ± -centers. It should be noted that, if = but t'^ ^ 0, the 
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eh-pair is locked in a two-site configuration. The two-electron(hole) hopping is governed 
by the transfer integral £jj that defines a probability amplitude for the exchange reaction: 

M ± + M T -> M T + M ± , 

or the motion of the electron (hole) center in the lattice of the hole (electron) centers. 
It is worth noting that in the conventional Hubbard-like models all the types of one- 
electron(hole) transport are governed by the same transfer integral: t'^ = t"j, while 
our model implies independent parameters for the disproportionation/recombination 
process and a simple quasiparticle motion on the lattice of M°-centers. In other words, 
we deal with a "correlated" single particle transport. 

2.2. Electron-hole Bose liquid, hard-core composite bosons and superconductivity in 
disproportionated systems 

Classical (or chemical) description of the mixed-valence systems implies a full neglect 
of the off- diagonal purely quantum CT effects: D^ 1 ' 2 ^ = t = 0, hence the valence of any 
site remains to be definite: 0, ±1, and we deal with a system of localized polar centers. 
Quantum description implies the taking account of the CT effects so that we arrive 
at quantum superpositions of different valence states resulting in an indefinite on-site 
valence and ionicity whose effective, or mean values (S z ) and (S 2 ) can vary from —1 to 
+1 and to +1, respectively. 

Simple uniform mean-field phases of the mixed-valent system include an insulating 
monovalent M°-phase (parent phase), mixed- valence binary (disproportionated) M ± - 
phase, and mixed- valence ternary ("under-disproportionated") M^-phase [I] . 

Insulating monovalent M°-phase with (S 2 ) = is a rather conventional ground state 
phase for various Mott-Hubbard insulators such as 3d oxides with a large enough positive 
magnitude of A > A cr parameter (U > 0). All the centers have the same bare M° valence 
state. Mixed-valence binary (disproportionated) M ± -phase with (S%) = 1 implies an 
overall disproportionation M° + M° — y M + M T . It is a rather unconventional phase 
for insulators. All the centers have the "ionized" valence state, one half the M + state, 
and another half the M~ one, though one may in common conceive of a deviation 
from the half-filling. A simplified "chemical" approach to M -phase as to a classical 
disproportionated phase is widely spread in solid state chemistry [13] . Typical ground 
state of such a classical phase corresponds to a checkerboard charge order, or longitudinal 
antiferromagnetic (staggered) pseudospin Ising-like Z-ordering. In systems with a short- 
range density-density coupling V nn the disproportionation phase transition M° — y M ± , 
or pseudospin reorientation occurs, if A < A cr = %V nn , where z is a number of nearest 
neighbours. In other words, the charge transfer instability of the parent M°-phase 
implies not only negative but also small enough positive A C t = Udd- 

The mixed valence M phase as a system of strongly correlated electron and hole 
centers appears to be equivalent to a Bose-liquid in contrast with the electron-hole 
Fermi-liquid in conventional semiconductors, hence it can be termed as electron-hole 
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Bose liquid (EHBL) [U \5\ dS]. Indeed, one may address the electron M~ center to be a 
system of a local composite boson (e 2 ) localized on the hole M + center: M~ = M + + e 2 . 

Three well known molecular-field uniform phase states of the M ± binary mixture, 
or EHBL phase can be specified as follows: 

i) charge ordered (CO) insulating state with (S z ) = ±1 and zero modulus of bosonic 
off-diagonal order parameter: ((S' 2 )! = 0; 

ii) Bose-superfluid (BS) superconducting state with (S z ) = 0, = e 2 ^; 

iii) mixed Bose-superfluid-charge ordering (BS+CO) superconducting state 
(supersolid) with < \(S g )\ < 1, (S 2 ) ^ 0. 

In addition, we should mention the high-temperature non-ordered (NO) Bose- 
metallic phase with <C S z ^>= 0. 

Disproportionation in such oxides as (Ca,Sr)Fe03, RMO3 yields an electron-hole 
Bose system with the 50% concentration of the composite bosons, or the half-filling. 
At variance with the parent quasi-2D cuprates the disproportionated 3d oxides from 
the very beginning can be addressed to be " negative- U" , or systems with A < A cr 
where the superconductivity can be driven by a deviation from the half-filling. 

How a typical insulating parent 3d M°-system, such as La2Cu04, can be driven to 
the disproportionated M , or EHBL phase? Simplest way is the hole/electron doping 
due to a nonisovalent substitution (NIS), that does solve two problems at once. First, 
it creates the impurity centers for a local condensation of the EH dimers, or pairs of 
M ± -centers thus shifting the system to the "negative-f/" regime. Second, the doping in 
cuprates such as La2_ a; Sr, r Cu04 and Nd2_a;Ce x Cu04 gradually shifts the M ± -phase away 
from half- filing making the concentration of the local S-bosons to be % = 0.5 — x/2 
(LSCO) or ub = 0.5 + x/2 (NCCO) thus promoting the superconductivity. In both 
the hole- and electron-doped cuprates we deal with composite S-bosons moving on the 
lattice of the hole centers CuOip, that makes the unconventional properties of the hole 
centers [29j EQ] to be common ones for the both types of cuprates. It is clear that the 
disproportionation scenario makes the doped cuprates the objects of bosonic physics. 
Obviously, the appearance of an inhomogeneous impurity potential under nonisovalent 
substitution deforms the phase diagram typical for uniform systems, in particular, in 
the " underdoped" regime where the system transforms from the parent phase M° into 
the disproportionated phase M . 

In the limit A — > — oo, the EHBL phase is equivalent to the lattice hard-core (he) 
Bose system with an inter-site repulsion whose Hamiltonian can be written in a standard 
form as follows (see Refs. [3~T| [32] and references therein): 

H hc = - E U 3 P{B\B 3 + B]B t )P + E VijNiNj - fi £ (3) 

i>j i>j i 

where P is the projection operator which removes double occupancy of any site. Here 
B\B) are the Pauli creation (annihilation) operators which are Bose-like commuting 
for different sites [B u B}] = 0, if i ^ j, [R h B\\ = \-2N u Ni = BjBf, N is a full number 
of sites. \x the chemical potential determined from the condition of fixed full number of 
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Figure 1. (Color online) The QMC T-n phase diagram of the 2D hard-core boson 
system (reconstruction of Fig. 2 from Ref. 34{, V\ = V nn ,V2 = V nnn )- Dotted 
line schematically circumscribes the region of the well developed BS+CO and BS 
superconducting fluctuations (local supercurrent response) giving rise to anomalous 
Nernst and local diamagnetism signals. 



N 

bosons Ni = J2 (Ni) or concentration n = Ni/N G [0, 1]. The t^ denotes an effective 

i=l 

transfer integral, Vij is an intersite interaction between the bosons. It is worth noting 
that near half-filling (n 1/2) one might introduce the renormalization N — > (A/j — 1/2), 
or neutralizing background, that immediately provides the particle-hole symmetry. 

The model of hard-core bosons with an intersite repulsion as a minimal model 
of the EH Bose liquid can be mapped to a system of pseudospins s = l/2 exposed to 
an external magnetic field in the Z-direction (Matsubara-Matsuda transformation [53]). 
For the system with a neutralizing background we arrive at an effective pseudo-spin 
Hamiltonian 

Hhc = Y, + + E J&m - A* E (4) 

i>j i>j i 

where J% = 2t y , = V lv s~ = j-B§+ = --$-B\ s z = -\+B\B u s± = T^(s x ±tP). 

In terms of the s — 1/2 pseudospins the non-ordered NO, or liquid phase, 
corresponds to a paramagnetic phase, Bose superfluid BS order does to a magnetic 
order in the XY plane, while charge density order CO does to a magnetic order in the 
Z direction. 

In Fig.[T] we present the phase diagram of the square lattice hc-boson Hubbard 
model with the nearest neighbour {nn) transfer integral t nn = t (the Josephson coupling) 
and repulsion V nn = 3t, derived from the quantum Monte-Carlo (QMC) calculations by 
Schmid et al. [M] (see, also Ref. [31]). Different filling points to CO phase, BS phase, and 
phase separated supersolid BS+CO phase. The AB line T KT (x) points to 2D Kosterlitz- 
Thouless phase transition; the C-B-D-C line points to the first order phase transition; 
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the D-E line T co (An) which can be termed as the pseudogap onset temperature T*(An) 
points to the second order Ising kind melting phase transition CO-NO into a nonordered, 
or normal fluid phase. At half-filling {tib = 0.5, An = 0) the system obviously prefers 
a checkerboard charge order (Neel antiferromagnetic order in the Z-direction below 
T co = 0.567 V nn given t = 0). 

It is worth noting that the phase diagram in Fig. [I] greatly resembles that of 
typical for doped quasi-2D cuprates [35] truly reproducing many important aspects of 
the normal and superconducting state, in particular, the pseudogap-like effects of the 
charge ordering and various signatures of the local superconductivity, that explain the 
anomalous Nernst [36] and local diamagnetism signals [37] . 

2.3. Topological phase separation in the 2D hard-core Bose system 

It is worth noting that all the lines in the phase diagram of the 2D hard-core Bose 
system in Fig. [T] point only to long range orders and do not concern the intricate intrinsic 
inhomogeneity. However, the CO state was shown [32] to be unstable with regard to 
a so-called CO+BS topological phase separation under doping, or deviation from the 
half-filling. The boson addition or removal in the half-filled hard-core boson system 
is assumed to be a driving force for a nucleation of a multi-center skyrmion-like self- 
organized collective mode that resembles a system of CO bubble domains with a Bose 
superfluid and extra bosons both confined in domain walls. The antiphase domain wall 
in the CO bubble domain appears to be a very efficient ring-shaped potential well for 
the localization of a single extra boson (or bosonic hole) thus forming a novel type 
of a topological defect with a neutral or single-charged (q=2e) domain wall. Such 
a topological CO+BS phase separation, more than likely being the dynamical one, 
rather than an uniform mixed CO+BS supersolid phase predicted by the mean-field 
approximation [31j, is believed to describe the evolution of hc-BH model away from 
half-filling. 

Fig. [2] presents the schematic view of the smallest neutral skyrmion-like bubble 
domain in a checkerboard CO phase for a 2D square lattice with an effective size of 
the three lattice spacings (~12A) [32]. The domain wall is believed to include as a 
minimum eight sites forming a ring-shaped system of four dimers each composed of two 
sites. There are two types of the domains which differ by a rotation by ±7r/2. Such 
a charged topological defect can be addressed as an extended (q = ± 2e) skyrmion-like 
mobile quasiparticle. On the other hand, the bubble seems to be a peculiar quantum of 
a "local superconductivity". 

The bubble domain may be addressed to be well isolated only if its close surrounding 
does not contain another domain(s) which could be involved to a "dangerous" frustrating 
overlap. Each domain in Fig]2]has z = 12 of such neighbors. Hence the concentration 
of well isolated domains can be written as follows: Pq(Au) = |An|(l — |An|) 2 , where z 
is a number of "dangerous" neighbors, An the doped boson concentration. The Po(An) 
maximum is reached at Ang = j^i • in our case ^ n o = 1/13, or « 0.077. With increasing 




Figure 2. Skyrmion-like bubble domain in the checkerboard CO phase with the 
eight site ring-shaped domain wall. Electron (M~) and hole (M + ) centers are marked 
by different filling. The four dimers within the domain wall are marked by I-IV. 
Their specific filling points to a bond-centered CO in the domain wall rather then a 
site-centered CO in the domain. Schematically shown are nearest and next-nearest 
neighbor domains which do not overlap with the central domain. Centers of all the 12 
nearest neighbor domains which could be involved to a " dangerous" frustrating overlap 
with the central domain are marked by crosses. 

doping the deviation of Po(An) from the linear law rises. On the other hand, knowing 
the effective domain area Sd ~ 9a 2 we can roughly estimate the limiting concentration 
of the single-domain model description to be An max « 1/9 ~ 0.11. 

The bubble's system is believed to reveal many properties typical for granular 
superconductors, CDW materials, Wigner crystals, and multi-skyrmion system akin in 
a quantum Hall ferromagnetic state of a 2D electron gas [32] • The concept of a static or 
dynamic bubble system introduces a large body of novelties into the physics of the EH 
Bose liquid related to a complex intra-bubble structure, bubble transport, bubble-bubble 
coupling, novel excitation modes, pinning of the bubble system etc. From a viewpoint 
of the novel charge order effects we point to an evolution of the bubble system upon 
lowering the temperature which implies a sequence of the isotropic liquid phase, the 
liquid- crystal phase, and the incommensurate bubble crystal phase with a quantum 
melting effect near the "magic" doping level, e.g. Arif, = 1/16, where the bubble 
lattice undergoes the structural phase transition. It is worth noting that the bubble 
crystallization is accompanied by different (pseudo)gap effects. From a viewpoint of 
the off-diagonal BS order we should point to a step-by-step formation of the local BS 
condensate starting with an intra-EH-dimer order at A;T~t nn , an intra-bubble order 



Perspectives of the disproportionation driven superconductivity 



12 



at fcT~t nnn , a short-range bubble order at kT~t bb (t 66 is a bubble transfer integral), 
formation of extended quasi- ID clusters which can undergo a transition to a phase- 
coherent state (filamentary or gossamer superconductivity [38j). Finally, upon the rise 
of the number of such clusters, the prerequisites for a percolative 2D and bulk 3D 
superconductivity are created. Attractively large temperatures of the emergence of 
the nonzero local BS condensate density engender different reasonable speculations as 
regards its practical realization. 

The topologically inhomogeneous phase of the hc-BH system away from the half- 
filling can exhibit the signatures both of s—,d—, and p- symmetry of the off-diagonal 
order [32] . The model allows us to study subtle microscopic details of the order 
parameter distribution including its symmetry in a real rather than momentum space, 
though the problem of the structure and stability of nanoscale domain configurations 
remains to be solved. 

It is worth noting that the bubble model predicts the optimal conditions for 2D and 
3D superconductivity in the doped cuprates at x opt ~2An ~ 0.154 and a suppression 
of superconductivity at x>2An max ^0.22. On the other hand, the bubble model 
predicts a so-called "self-doping" (SD) mechanism of superconductivity in "negative- 
U" , or M ± systems. Such a mechanism is realized when the CO phase in the system 
appears to be unstable with regard to a nucleation of positively and negatively charged 
bubbles preserving a charge neutrality. Probably, such a "self-doping" explains the 
superconductivity in YBa 2 Cu 4 8 [39] and LiFeAs [4"0] . 

It should be noted that the bubble system can manifest two types of the dispersive 
features that reflect a dominance of the kinetic energy (the bubble transfer) or the 
potential energy (the bubble-bubble coupling), respectively. Indeed, a single mobile 
bubble has energy e(q), while different static bubble crystal configurations (Wigner 
electronic crystals) can be specified by the quasimomentum q and corresponding energy 
E(q). 

3. Disproportionation and local composite bosons in 3d 71 systems 

3.1. Local composite bosons in 3d 11 systems 
As we have seen above the disproportionation 

3d n + 3d n -> 3d n+1 + 3d n ~ l 

is addressed to be a candidate mechanism to drive the superconductivity in nominally 
insulating 3d compounds. A simple view of the disproportionated system implies 
the electron(hole) center 3d n+l (3d n ~ 1 ) to be composed of hole (electron) center plus 
electron(hole) coupled pair 3oP(3d 2 ) that seems to constitute a composite electron(hole) 
local boson. In other words, the disproportionated system is anticipated to be a 
system of local bosons moving in a lattice formed by electron(hole) centers. The 
disproportionation scenario implies the composite local boson to be anyhow an integral 
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part of a stable many-electron atomic configuration, 3<i n+1 or 3<i n_1 . Microscopic nature 
of the "attractive force which could overcome the natural Coulomb repulsion between 
two electrons which constitute a Bose pair" [H] evolves both from main intra-atomic 
correlations and electron-lattice polarisation effects [22] providing the stability of 3d n±l 
configurations. 

However, for the representation of local composite bosons to justify we should take 
into account some rather strict limitations imposed by the orbital structure of the 3d 
states and the genealogy of many-electron 3d n,n±l states. First of all, the effective 
two-particle two-site transfer integrals 

(3dJ +1 3dr 1 |#|3dT 1 3d2 +1 ) > 
that describe the exchange reaction 

3d n+1 + 3d" 1 - 1 -)• 3d"- 1 + 3d n+1 

are due to reduce to a simple form 

(3d? +1 3dr 1 |#|3dr 1 3d£ +1 ) = (3dl\H\3d 2 2 ) = t 12 

we need to introduce directly the local boson representation. The two-particle transfer 
integral t±2 is a main parameter that governs both the transport and superconducting 
perspectives of the disproportionated system. It can be written as follows: 

11 A-dd 

where the first term describes a simultaneous tunnel transfer of the electron pair due 
to the Coulomb coupling V ee and may be called as a "potential" contribution, whereas 
the second describes a two-step (20-11-02) electron-pair transfer via successive one- 
electron transfer due to the one-electron Hamiltonian h, and may be called as a "kinetic" 
contribution. The value of the composite boson transfer integral is closely related to 
that of the exchange integral, i.e. t i2 ~ Jyi- 

Justification of the local composite boson implies a specific genealogy of the many- 
electron states ^(3rf n+1 ) and ^(3rf n_1 ): the wave functions should have a simple 
structure 

q(3d n+1 ) = q(3d 2 )^(3d n - 1 ), ^^d 71 - 1 ) = ^(3d 2 )^(3d n+1 ) (5) 

in electron, or hole representation, respectively. To gain more insight into the main 
features of the disproportionation reaction hereafter we analyze the d-d CT transition 
3d n + 3d n — y 3d n+l + 3d n ~ l in strongly correlated 3d systems with a high (octahedral, 
tetrahedral, and cubic) local symmetry. For crystal field configurations such as ^e™ 2 the 
optimal disproportionation scheme (jSJ) implies single configurations both for \I/(3<i n+1 ) 
and ^{3d n ~ l ). It is easy to see that the optimal conditions for the superconductivity 
are expected for parent 3d n systems with 3<i n±1 configurations which correspond either 
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empty, or filled and half-filled t 2g and e g orbitals. Second, to minimize the reduction 
effect of the electron-lattice interaction and avoid the localization, we need the 5-type 
(Ai g ,A2g) orbital symmetry of the local boson that provides the conservation of the 
orbital degeneracy when it moves on the lattice. Third, to minimize the reduction effect 
of the spin degrees of freedom, we need spin-singlet local bosons. Fourth, to maximize 
the transfer integral we need a participation of the strongest a bonds. In other words, 
we should point to a specific role of the one-electron orbitals forming the two-electron 
boson state. Furthermore, the optimized ground state of the local boson should be 
well isolated, that is separated from excited states by a gap on the order of 1 eV. In 
any case, along with the electron-lattice polarisation effects which stabilize the 3d n±l 
configurations and minimize the CT energy, the many-electron configuration, intra- 
atomic correlations, force and symmetry of the crystal field seem to be most important 
factors of the disproportionation scenario. All the above constraints restrict strongly 
the number of 3d n systems to be candidates for the high-T c superconducting materials 
in frames of the disproportionation scenario. 

3.2. Optimal disproportionation schemes for different 3d 11 ions in high- symmetry 
crystal field 

Below we address the optimal disproportionation schemes for different 3d n ions in high- 
symmetry crystal field. It is worth noting that for octahedral or tetrahedral (cubic) 
crystal field the electron filling starts with the t 2g or e g orbitals, respectively. Depending 
on the relation between the crystal field and intra-atomic correlation energies we will 
consider either high-spin (HS) or low-spin (LS) complexes such as MeOg or MeC^ in 3d 
oxides. 

3d 1 ions (Ti 3+ , V 4+ ). Disproportionation 3d 1 +3d 1 — > 3d°+3d 2 gives rise to a 
system of local composite bosons, or electron pairs 3d 2 , moving in a lattice formed by 
the 3d centers. Simplest electronic configuration 3d 1 is realized as t\ g configuration 
for octahedral or e g configuration for the cubic or tetrahedral crystal field. For the t\ g 
configuration we arrive at an orbital and spin triplet ground state 3 T lg with seemingly 
no optimistic expectations as regards the superconductivity. Some perspectives may 
be related to a spin-orbital coupling Vso that stabilizes a singlet state with a net 
momentum J=0, though the electron-lattice interaction can destroy spin-orbital coupling 
giving rise to an Jahn- Teller polaron. More attractive expectations are related to the 
e 2 configuration in the tetrahedral crystal field when we arrive at an S-like though spin- 
triplet ground state 3 A 2g . In other words, there appears an opportunity to realize an 
unconventional system of spin-triplet local S-like bosons moving on the electronically 
"inactive" lattice (see Fig. [I]). 

3d 2 ions (V 3+ , Cr 4+ ). Two-electron configuration 3d 2 is realized as t\ g 
configuration for octahedral or e 2 configuration for the cubic or tetrahedral crystal field. 
Disproportionation 3d 2 +3d 2 — > 3d x +3d 3 seems to gives rise to a system of electron pairs 
3d 2 , or local composite bosons, moving in a lattice formed by the 3d 1 centers. However, 
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the genealogy of the low-spin e 3 ; 2 E g configuration does not permit the insertion of a well 
defined composite boson. Indeed, the genealogical link of the e 3 and e l g configurations 
reads as follows: 



e 3 - 2 E 

c gi ^9 



1 

72 



^\{e^A lg -e]-. 2 E g )-±= (e^E n ; el : 2 E, 



J gl ^gi °# • ^9 



elYA 2g ;el: 2 E g ), (6) 



and implies a participation of not one but three different terms of the e configuration. 



9 

Genealogy of the high-spin t\ g ; A A 2g configuration obeys the master condition (see 
Exp. (2)) that permits the formation of a well defined composite boson with HS- 
configuration 3 T lg . It means that the disproportionation reaction 

±2 _|_ ±2 . rl , f3 
L 2g ' l 2g ^ L 2g ' L 2g 

gives rise to a system of the local electron type bosons (t 2 g ) on the lattice with the on-site 
t\ g configuration or hole type bosons (t^) on the lattice with the on-site HS-t| 9 ; A A 2g 
configuration. However, in any case such a composite boson appears to be charged by 
both spin and orbital degrees of freedom, in addition, this is prone to move on the lattice 
with a spin or spin and orbital degeneracy that strongly restricts "superconducting 
perspectives" . Some hopes may be related to a spin-orbital coupling Vso that stabilizes 
a singlet state with a net momentum J=0, though the electron-lattice interaction will 
compete with the spin-orbital coupling giving rise to an Jahn- Teller polaron. 

3d 3 ions (V 2+ , Cr 3+ , Mn 4+ ). Three-electron configuration 3d 3 is realized as t\ g 
configuration for octahedral, e 3 or e 2 t 2g configuration for strong or weak tetrahedral 
or cubic crystal field, respectively. However, the local boson representation can be 
introduced only for tetrahedral/cubic high-spin configuration: 

e 2 g t 2g + e 2 g t 2g ^ e 2 + e 2 g tj g 



when the disproportionated system can be addressed to be a system of t 2 g ; 3 T l9 local 
composite bosons on the lattice with the on-site e 2 ; 3 A 2g 3d-configuration. Obviously, 
we arrive at delusive perspectives for the high-T c superconductivity. 

3d 4 ions (Mn 3+ , Cr 2+ , Fe 4+ ). The electron configuration near the half-filling, 
3d 4 , is particularly interesting given the octahedral crystal field. Then the high-spin 
(Hund) 3d A configuration disproportionates as follows: 

{t%e\- 5 E g ) + {tl g e\- 5 E g ) -> {t\ g] 4 A 2g ) + (4 g e 2 g] 6 A lg ) , (7) 

f 3 2. 6 a \ _ (+3 . 4 A \ v ( a 2\. 3 



where (tl g e 2 ,Ai g ) = (t 2g ; A 2g ) x (e 2 ); A 2g . In other words, the disproportionated 
system can be viewed as a system of spin-triplet composite electron bosons (e 2 ); 3 A 2g 
moving in the lattice composed of the localized hole spin-3/2 S-type {t\ g ; 4 A 2g ) centers 
(see Fig.[T]). Interestingly that we can address the system in other way to be a system 
of spin-triplet composite hole bosons (e^); 3 A 2g moving in the lattice composed of the 
localized electron spin-5/2 S-type (t^e 2 ; 6 Ai g ) centers. It is worth noting that the 
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electron-type boson obeys a conventional (ferromagnetic) Hund rule on lattice sites while 
the hole-type boson obeys an unconventional (antiferromagnetic) "anti"-Hund rule on 
lattice sites. Local composite bosons can be formed also for the low-spin octahedral 
or high-spin tetrahedral d 4 systems, however, we arrive at the composite boson with a 
3 Ti g or t 2 ; 3 ^~ig configuration and a spin-3/2 or 5/2 lattice. 

3d 5 ions (Fe 3+ , Mn 2+ ). The low-spin 3d 5 configuration appears to be near filling 
both in octahedral and tetrahedral crystal field when it disproportionates as follows: 

(4) + (4) ^(4;%) + (4;^) (s) 

or 

(44) + (44) (ej; %,) + (44; 3 T lg ) (9) 

where {t\ g) 3 T lg ) = 3 T l9 ) x (t% g ] l Ai g ). In other words, the disproportionated system 
given the octahedral crystal field can be viewed as a system of spin and orbital triplet 
composite hole bosons (i| fl ; 3 T lg ) moving in a lattice composed of the spinless S-type 
((4> l Ai g ) centers. Given tetrahedral crystal field we arrive at a system of spin and 
orbital triplet composite electron bosons {t% g ; 3 Ti g ) moving in a lattice composed of the 
spinless S-type ((e 4 ; 1 A lg ) centers. Orbital degeneracy of the local bosons can give rise 
to the formation of the JT polarons followed by its localization. 

3d 6 ions (Fe 2+ , Co 3+ ). In this case the specific genealogy flSJ) of polar 3d n±1 
centers can be realized only for high-spin (Hund) configuration. In tetrahedral or cubic 
crystal field the 3d 6 system disproportionates as follows: 

(44; + (44; (44; + (44; 4 ^ 9 ) , (10) 

where (e 2 g t\ g -^ Ai g ) = (e 4 t 3 , g ; 4 A2 9 ) x (e 2 ); 3 v4 2 c,. In other words, the disproportionated 
system can be viewed as a system of spin-triplet composite hole bosons (e 2 ); 3 A 2g moving 
in the lattice composed of the electron S=3/2 (e 4 t\ g) 4 Ai g ) centers. We can address the 
system in other way to be a system of spin-triplet composite electron bosons (e 2 ); 3 A 2g 
moving in the lattice composed of the localized hole spin-5/2 S-type (44' cen t ers 
(see Fig.[T]). The situation resembles that of the Hund 3c? 4 electron configuration given 
the octahedral crystal field. 

In octahedral crystal field the HS-3g! 6 system disproportionates as follows: 

(44; + (44; -> (44; + (44; %) • ( n ) 

The disproportionated system can be viewed as a system of spin-triplet composite 
electron-type bosons (t\ g ; 3 Ti g ) moving in the lattice composed of the hole S=5/2 
t^e 2 ; 6 A lg centers. Orbital degeneracy for the composite boson implies a competition 
of the vibronic JT coupling and the kinetic energy. 

3d 7 ions (Ni 3+ , Co 2+ ). For low-spin (non-Hund) configuration in octahedral 
crystal field the 3d 7 system disproportionates as follows: 

(44; 2 E g ) + (44; 2 E g ) -> (t%; % g ) + (44; 3 A 2g ) , (12) 
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Figure 3. (Color online) Disproportionation schemes for the high-symmetry crystal 
held which are believed to provide optimal conditions for the disproportionation driven 
superconductivity. 



where it\ g e 2 g \ 3 A 2g ) = {t% g \ 1 A\ g ) x (e 2 )] 3 A 2g . In other words, the disproportionated 
system can be viewed as a system of spin-triplet composite electron bosons (e 2 ); 3 A 2g 
moving in the lattice composed of the hole S=0 (t 2g ', 1 Ai g ) centers (see Fig. [I]). 

For high-spin (Hund) configuration in octahedral crystal field the 3d 7 system 
disproportionates as follows: 

(44; 4T ^) + (44; %) -> (44; + (44; 3 ^.) , (is) 

The disproportionated system can be viewed as a system of spin-triplet composite hole 
bosons t% g ] 3 Ti g moving on the lattice composed of the electron S=l (44' 3j ^2 g ) centers. 

3d 8 ions (Ni 2+ , Cu 3+ ). For the ground state configuration in octahedral crystal 
field the 3d 8 system disproportionates as follows: 

{t\/ g - 3 A 2g ) + (f 2g e 2 g - 3 A 2g ) -> (4eJ; 2 E g ) + (t%e 3 g , 2 E g ) , (14) 

with the formation of the two Jahn- Teller centers which are prone to a vibronic 
localization. However, the genealogy of the e 3 ; 2 E g configuration does not permit the 
insertion of the well defined composite boson (see Exp.[H]). The 3d 8 system given the 
tetrahedral crystal field disproportionates as follows: 

(44; 3T hv) + (44; (44; 4 ^,) + (44; 2T ^) > ( 15 ) 

In other words, such a disproportionated system can be viewed as a system of spin- 
triplet composite electron- type JT bosons [t 2g ; 3 Ti g ) moving on the lattice composed of 
the hole spin-3/2 (e^| g ; 4 A 2g ) S-type centers. 

3d 9 ions (Cu 2+ ). In octahedral crystal field the 3c? 9 system disproportionates as 
follows: 

(t%e 3 g ; 2 E g ) + (t b 2g e 3 g ; 2 E g ) -> {t%e 2 g] 3 A 2g ) + {t%e g] l EA lg ) , (16) 
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where (4 e <?' 3 ^2 9 ) = (^2 9 e <j> 1 Ai g ) x (e 2 ,); 3 ^^- in other words, the disproportionated 
system can be viewed as a system of spin-triplet composite hole bosons (e 2 ); 3 A 2g moving 
in the lattice composed of the electron spin- (t| fl ej; 1 A lg ) centers (see Fig.[TJ). 

In tetrahedral or cu&ic crystal field the 3d 9 system disproportionates as follows: 

(44; 2 T 2g ) + (e%; 2 T 2g ) (ejtj,; 3 T l9 ) + , (17) 

where (e 4 ,^; 3 ^is) = ( e ^2 g ! x GSjg)) 3 ^- I 11 other words, the disproportionated 

system can be viewed as a system of spin-triplet composite hole bosons 3 Ti g moving 
in a lattice composed of the electron spinless {e^t 2g ; centers. 

5.5. Summary 

Let's summarize our findings. 

i) The most effective e 2 g configuration of the local composite boson providing 
maximal values of the boson transfer integral without strong reduction effects of the 
electron-lattice coupling is realized only for several optimal d n configurations. These 
are HS-c? 4 , LS-d 7 , d 9 configurations given octahedral crystal field, and d 1 , RS-d 6 
configurations given tetrahedral crystal field |42j ; 

ii) All these bare systems are characterized by E-type orbital degeneracy, i.e. these 
are prone to a strong Jahn- Teller effect. In all the instances the disproportionation 
reaction lifts the bare orbital degeneracy, that is it has a peculiar "anti- Jahn- Teller" 
character [33] . 

iii) All these bare systems are characterized by a strong suppression due to a vibronic 
reduction of the one-particle (electron or hole) transport, while after disproportionation 
these are characterized by an effective two-particle (local boson) transport. 

iv) Undesirable spin-triplet structure of the local composite bosons appears to be 
a common feature of all the above mentioned disproportionated systems with a high- 
symmetry crystal field. 

v) At variance with d l , d 7 , and d 9 systems the d A and d 6 systems reveal unavoidable 
coexistence of the spin-triplet bosons and a magnetic lattice. 

All these results are summarized in Table I where we have added examples of 3d 
systems with optimal 3d n configurations. The Fig. [3] illustrates the final electronic 
configurations and two-particle e 2 g transport for the disproportionation reaction in 
optimal <i 4 ' 6 ' 7 ' 9 systems. 

Our consideration does support and generalize a purely electronic motivation for 
"negative U" centers which was based upon a stabilization of the electron-rich closed- 
shell configurations such as ns 2 or nd 10 through the charge excitation effects or of 
half-filled configurations such as nd 5 through the exchange-correlation effects (Hund's 
rule) The approach does explain the missing oxidation states of ions in s 1 electronic 
configuration and predict the disproportionational, or charge transfer instability in 
systems with ions having nominally ns 1 (Hg + , Tl 2+ , Pb 3+ , Bi 4+ ) , d 9 (Cu 2+ ), d A (Mn 3+ , 
Fe 4+ ), d 6 (Fe 2+ , Co 3+ ) valence electronic configuration. 
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Table 1. 3d™ JT-systems optimal for the disproportionation driven superconductivity 
(see text for detail) 
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In the above we proceed with a conventional approach to the electronic structure 
of the 3d compounds that implies a predominantly 3d antibonding character of the 
molecular t 2 g and e g orbitals in octahedral M06 or tetrahedral MO4 clusters. Electron 
and spin density in these clusters is distributed between the 3d ion and ligands 
that should be taken into account, in particular, when address a so-called "partial" 
disproportionation such as Ni^ 3+5 ^ + - Ni^ 6 ^ in RNiC>3[12j. However, upon lowering 
the crystal field symmetry and anomalous strengthening of the tip -hybridization in the 
end of the d n series the situation can change in several important points. Remarkable 
example is provided by quasi-2D cuprates where we proceed with a strong axial 
distortion of the Cu06 octahedra, actually with square CUO4 plaquettes (D^h point 
symmetry) as a basic element of the crystal and electronic structure. Single hole 
&i 5 (oc <i :r 2_j / 2)-state of Cu 2+ ion is a typical one for the square-planar coordination 
of oxygen ions whereas two-hole Cu 3+ state with the same coordination of oxygen 
ions seldom exists. Instability of Cu 3+ ions with 3d 8 , or two-hole 3d 2 configuration 
is obviously related to a strong inter-electron(hole) repulsion. However, the two-hole 
state in CuO^ - center to be a cluster analogue of Cu 3+ ion can be stabilized by a 
location of the additional hole to low-energy predominantly oxygen O 2p molecular 
orbitals thus providing a sharp suppression of the inter-electron(hole) repulsion at a 
relatively small loss in the one-hole energy. In 1988 Zhang and Rice jH] have proposed 
that the doped hole forms a well isolated local spin and orbital 1 Ai g singlet state which 
involves a phase coherent combination of the 2p<r orbitals of the four nearest neighbor 
oxygens with the same h\ g symmetry as for a bare Cu 3d x 2_ y 2 hole. The spin and 
orbital singlet ]? Xg \ 1 A lg , or Zhang-Rice singlet, may be viewed as a local hole composite 
boson. In other words, the disproportionation 3d 9 + 3d 9 — > 3d 10 + 3d 8 in the tetragonally 
distorted octahedral 3d 9 = 3d 1 (Cu 2+ ) systems gives rise to a system of composite hole 
S=0 bosons whose spin and orbital structure provides maximal values of the effective 
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boson transfer integral and minimal boson effective mass, that is optimal conditions for 
the disproportionation driven high-T c superconductivity. Such a situation is unlikely 
realized for the d l configuration in a tetrahedral crystal field because of a weak p- 
d covalency. Thus, our analysis shows that in fact only d 9 system can be a major 
candidate for maximal T c 's. 

Above we did not concern the energy stability of the final disproportionated 
"negative-?/" state. It should be noted that together with a sizeable recombination 
energy this implies well isolated ground states of d n±1 configurations. 

4. Electron-hole Bose liquid with spin-triplet s = 1 local composite bosons 

Above we have presented a unified approach to the disproportionation phenomenon in 
different 3d compounds though we did not touch upon such important points as electron 
and lattice polarization effects, vibronic coupling, and the role played by the spin and 
orbital degrees of freedom. All these can strongly deform the phase diagram of the 3d 
compounds from the predictions of the purely charge order parameter model. 

Minimal model of the EHBL phase does not imply intervention of the spin degrees 
of freedom only in quasi-2D cuprates. Indeed, the model considers such a cuprate 
to be a system of the the spin and orbital singlet l Ai g local S-bosons moving on the 
lattice formed by hole centers with the well isolated spin and orbital singlet Zhang- Rice 
l A\ g ground state. For all other "optimal" 3d n systems listed in Table I we arrive at an 
unavoidable spin-triplet s — 1 structure of the composite local bosons moving on the spin 
(bare octa-HS d 4 and tetra-HS d 6 configurations) or spinless (bare tetra- d 1 and octa- 
LS d 7 configurations) lattice. In the absence of the external magnetic field the effective 
Hamiltonian of such an electron-hole Bose liquid takes the form of the Hamiltonian of 
the quantum lattice Bose gas of the triplet bosons with an exchange coupling [19] : 

H = Hqlbg + H ex = ^2 t B(ii)B\ m Bj m 
+ £ V tjni n, - ^ £ ni + E J ?( S * ■ Si) 

i>j i i>j 

+ E -Si)+E J$(* ■ ^) + E ■ ^ ■ ( 18 ) 

Here B\ m denotes the s = 1 boson creation operator with a spin projection m at the 
site i] Bi m is the corresponding annihilation operator. The boson number operator 
him = B\ m B im at z-site due to the condition of the on-site infinitely large repulsion 
Vu +oo [hardcore boson) can take values or 1. 

The first term in ( fl8l) corresponds to the kinetic energy of the bosons, tsiij) is 
the transfer integral. The second one reflects the effective repulsion (V^ > 0) of the 
bosons located on the neighboring sites. The chemical potential \x is introduced to fix 
the boson concentration: n = Y^i(fk) ■ F° r EHBL phase in parent system we arrive at 
the same number of electron and hole centers, that is to n = ^. The remaining terms in 
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(I18p represent the Heisenberg exchange interaction between the spins of the hole centers 
(term with J hh ), spins of the hole centers and the neighbor boson spins (term with J hb ), 
boson spins (term with J bb ), and the very last term in ( Tl8l) stands for the intra-center 
Hund exchange between the boson spin and the spin of the hole center. In order to 
account for the Hund rule one should consider J^ b to be infinitely large ferromagnetic. 

Generally speaking, this model Hamiltonian describes the system that can be 
considered as a Bose-analogue of the one orbital double-exchange model system [35] . 

More or less simple manifestation of spin degrees of freedom are anticipated for 
tetra-c? 1 and LS-octa-c? 7 configurations where we proceed with spin-triplet composite 
bosons, however, moving on the spinless lattice. 

Estimates for different superexchange couplings given the bond geometry typical 
for manganites such as LaMnO"3 [19] predict antiferromagnetic coupling of the nn hole 
centers (J hh > 0), antiferromagnetic coupling of the two nearest neighbor bosons 
(J bb > 0), and ferromagnetic coupling of the boson and the nearest neighbor hole 
centers (J hb < 0). In other words, we arrive at highly frustrated system of triplet 
bosons moving in a lattice formed by hole centers when the hole centers tend to order 
G-type antiferromagnetically, the triplet bosons tend to order ferromagnetically both 
with respect to its own site and its nearest neighbors. Furthermore, nearest neighboring 
bosons strongly prefer an antiferromagnetic ordering. Lastly, the boson transport prefers 
an overall ferromagnetic ordering. 

Competition of the Heisenberg exchange and the bosonic double exchange in 
"optimal" systems such as ferrates (Ca,Sr)Fe03 or manganites RMnOs with bare octa- 
HS d 4 configurations of the transition metals can be easily demonstrated on an example 
of the d 5 -d 3 pair, or EH-dimer. The net spin of the EH-dimer is S = Si + S2, where 
S x (Si = 5/2) and S 2 (S 2 = 3/2) are spins of Fe 3+ and Fe 5+ (Mn 2+ and Mn 4+ ) ions, 
respectively. In nonrelativistic approximation the spin structure of the EH-dimer will 
be determined by isotropic Heisenberg exchange coupling 

V ex = J(Si • S 2 ), (19) 

with J being an exchange integral, and the two-particle charge transfer characterized 
by a respective transfer integral which depends on spin states as follows: 

(~; SM\H B \ 3 -^, SM) = ±S(S + l)t B , (20) 

where t B is a spinless transfer integral. Making use of this relation we can introduce an 
effective spin-operator form for the boson transfer as follows: 
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2(S 1 -S 2 ) + S 1 (S 1 + l) + S 2 (S 2 + l) , (21) 



which can be a very instructive tool both for qualitative and quantitative analysis 
of boson transfer effects, in particular, the temperature effects. Both conventional 
Heisenberg exchange coupling and unconventional two-particle bosonic transfer, or 
bosonic double exchange can be easily diagonalized in the net spin S representation 



Perspectives of the disproportionation driven superconductivity 



22 




Figure 4. (Color online) Spin structure of the EH-dimer, or self-trapped CT exciton 
with a step-by-step inclusion of one- and two-particle charge transfer. Arrows point to 
electric dipole moment for bare site-centered dimer configurations. 



so that for the energy we arrive at 

Es = ±[S(S + 1) - y] ± ^S(S +l)t B , (22) 

where ± corresponds to two quantum superpositions |±) written in a spin representation 
as follows 

\SM) ± =-±=(\^SM)±\^SM)), (23) 

with s- and p-type symmetry, respectively. It is worth noting that the bosonic double 
exchange contribution formally corresponds to a ferromagnetic exchange coupling with 
Jb = — iq|^b|- 

We see that the cumulative effect of the Heisenberg exchange and the bosonic double 
exchange results in a stabilization of the S = 4 high-spin (ferromagnetic) state of the 
EH-dimer provided |te| > 10 J and the S = l low-spin (ferrimagnetic) state otherwise 
(see Fig.H]). As in conventional (one-electron) double exchange model [16] the latter 
frustration between the Heisenberg exchange and the bosonic double exchange stabilizes 
a noncollinear ordering in bulk systems. 

5. Discussion 

Our approach provides a comprehensive understanding of the well established 
disproportionation in ferrates (Ca,Sr)Fe03 and nickellates RNiOs with octa-HS d A and 
octa-LS d 7 configurations of the transition metals. It is worth noting that the flexible 
perovskite structure of both compounds facilitates a screening of the effective correlation 
parameter Udd- In (Ca,Sr)Fe0 3 the effect is believed to strengthen due to a large enough 
polar izability of the Ca 2+ or Sr 2+ ions while in RMO3 we deal with a sizeable reduction 
of the effective Udd due to a strong Ni-0 covalency. 

An incommensurate helicoidal spin ordering observed both in CaFeOs and SrFeC>3 
up to very low temperatures can be explained as a result of a competition between 
conventional exchange coupling and the bosonic double exchange in EH dimers 
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M ± -H-M =F . Interestingly, the helicoidal spin ordering gives rise to a full (SrFeOs) or 
partial (CaFeOa) suppression of the charge order resulting in a metallic or insulating 
temperature dependence of the resistivity, respectively jH]. In any case both theoretical 
and experimental study of the phase diagram for (Ca,Sr)Fe0 3 system deserves further 
work, especially, aimed to a search of a possible superconductivity 

Nominal electronic HS-octa-configuration d A of the Fe 4+ ion is analogous 
to that of Mn 3+ ion in perovskite (ortho)manganites RMn03 that implies 
similar disproportionation features. Indeed, similarly to ferrates (Ca,Sr)Fe03 the 
orthomanganites RMnC>3 reveal a high-temperature metallic phase which may be 
associated with a fully disproportionated R(Mn 2+ Mn 4+ )C>3 phase [T8l [19] . or EHBL 
system. Such a conclusion is particularly supported by an anomalously small magnitude 
of the thermopower. Interestingly that under spin ordering this manganite phase in 
contrast with the ferrates could become a ferromagnetic EHBL metal as bosonic double 
exchange in manganese EH dimers overcomes the Heisenberg exchange in favour of 
the ferromagnetic spin state [IH] . In other words, in manganites we deal with the 
high-spin ferromagnetic ground state of the EH dimer (see Fig.Hk), while in ferrates 
with the low-spin ferrimagnetic one (see Fig. lib). However, actually upon lowering 
the temperature one observes a first order phase transition at T=T jt {Tjt ~ 750 K in 
LaMnOs) from the high-temperature fully disproportionated metallic EHBL phase to 
a low-temperature orbitally ordered insulating phase with a cooperative Jahn- Teller 
ordering of the occupied e g -orbitals of the Mn 3+ 6 octahedra accompanied by A- 
type antiferromagnetic ordering below (Tjv~140K in LaMnOa) [T8| IT9] . In other 
words, the vibronic (Jahn- Teller) coupling in bare (parent) Mn 3+ phase suppresses the 
anticipated charge order and does promote a recombination M — > M° transition. 

Interestingly that the nonisovalent substitution and/or nonstoichiometry seems 
to revive the disproportionated phase and such manganites along with a metallic 
ferromagnetism with a colossal magnetoresistance reveal many properties typical for 
superconducting materials. These are a steep jump of the specific heat, gap opening 
in tunneling, and giant proximity effect. Several manganite samples were shown to 
exhibit a negative diamagnetic susceptibility in a wide temperature range up to room 
temperature [48 j . To explain many experimental manifestations of superconductivity 
Kim [19] has proposed a superconducting condensation in manganites. In this 
scenario, colossal magnetoresistance in Lai_ x Sr x Mn03 is naturally explained by the 
superconducting fluctuations with increasing magnetic fields. This idea is closely related 
to the observation of anomalous proximity effect between superconducting YBCO and 
a manganese oxide, Lai-^Ca^MnOs or Lai_ z Sr x Mn03 [50], and also the concept of a 
local superconductivity manifested by doped manganites [51] . 

A strong evidence of the revived EHBL phase below first order phase transition 
M ± — y M° from the high-temperature EHBL phase to low-temperature AFM insulating 
phase is obtained very recently by Nath et al. [52] . The authors observed an electric field 
driven destabilization of the insulating state in nominally pure LaMnC>3 single crystal 
with a moderate field which leads to a resistive state transition below 300 K. The effect 
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has been explained as bias driven percolation type transition between two coexisting 
phases, where the majority phase is a charge and orbitally ordered polaronic AFM 
insulating phase and the minority phase is a bad metallic EHBL phase |52j . 

Very recently distinct signatures of high-temperature disproportionated EHBL 
phases are revealed also in other manganites, such as LaMn70i2 [53] with quadruple 
perovskite structure and YBaMn 2 06 [54J. A dramatic suppression of the Curie- Weiss 
susceptibility in the undistorted phase of LaMn 7 0i2 above Tjt ~ 650 K is explained 
by a picture of long-lived low-spin S=l EH dimers (see Fig.Hb) formed by the charge 
transfer between neighboring Mn sites [53] . High-temperature electron spin resonance 
studies of A-site ordered mixed- valent manganite YBaM^Oe [HI] have revealed two 
puzzles. First, in the metallic high-temperature phase the ESR signal stems only from 
all Mn 4+ core spins. Second, the ESR linewidth obeys the Korringa-like linear high- 
temperature increase over a wide temperature range of 400 K clearly pointing to a truly 
metallic state. Both findings strongly support our suggestion of the high-temperature 
EHBL phase of manganites to be a system of spin-triplet composite e 2 g bosons moving 
on the lattice of spin-3/2 Mn 4+ centers. Furthermore, one may speculate that the 
Korringa-like behavior of the ESR linewidth could be a signature of such a charge- 
disproportionated picture [51] . Our prediction of the disproportionation driven high-T c 
superconductivity in the HS-c/ 4 octa-systems needs in a more direct validation for the 
hole doped manganites such as Lai^Sr^MnOs. 

High-temperature metallic phase of rare earth nickelates RNIO3 with LS-octa-(i 7 
configuration of Ni m+ ions resembles that of orthomanganites, however, with a Pauli- 
like magnetic susceptibility and evidences a full d-d disproportionation. Yet the low- 
temperature behavior of manganites and nickellates strongly differs. For rhombohedral 
LaNiC>3 the system remains metallic down to 1.5 K, whereas orthorhombic RMO3 (R 
=Pr,... Lu) exhibit a first order phase transition to a charge ordered insulating state 
upon cooling below Tco spanning from 130 K for Pr to ~ 550-600 K for heavy rare earths. 
Both the lighter and heavier orthorhombic RMO3 compounds are very similar from the 
point of view of their local electronic and magnetic CO state despite the strong change of 
the metal-to-insulator NO-CO transition temperature. All these exhibit clear signatures 
of the charge disproportionated state with two types of Ni centers corresponding to 
alternating large NiOg 0- (Ni 2+ center) and small NiOg - (Ni 4+ center) octahedra strongly 
differing in magnetic moments (~2/z# and ~0, respectively). Regarding the magnetic 
order, these nickellates adopt an antiferromagnetic structure defined by the propagation 
vector (1/2,0,1/2) at a transition temperature TV spanning from 130 K to ~200K. 
Antiferromagnetic ordering can be explained by a rather strong superexchange nnn 
(next-nearest neighbor) coupling of magnetic S=l Ni 2+ centers. Unfortunately, the 
antiferromagnetism stabilizes the charge order and suppresses the superconductivity. 

At variance with the charge transfer unstable perovskites ((Ca,Sr)Fe03, RMnOs, 
RMO3) the quasi-2D nickellates AMO2 (A = Ag, Li, Na) reveal existence of 
unconventional ground states stabilized by the frustrated triangular lattice geometry 
from a cooperative JT ordering of Ni 3+ ions in NaNi02 to a moderately charge ordering 
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3Ni m+ — > Ni 2+ +2 Ni 3 5+ in antiferromagnetic metal AgNi0 2 [53]. In the case of LiNi0 2 
there could be a competition between charge ordering and orbital ordering for the ground 
state, the nickel valency could be a mixture of 2+, 3+, and 4+ [56J. 

In all the CT-unstable 3d compounds with optimal electron configurations 
addressed above the JT coupling in nearly octahedral complexes (FeOg, MnC>6, NiOe) 
competed with the charge disproportionation and in several compounds this determined 
the orbitally ordered ground state and the low-temperature phase thus suppressing 
the anticipated superconductivity. The JT coupling reaches its maximal magnitude in 
Cu 2+ octahedral complexes that seemingly should inhibit the disproportionation in Cu 2+ 
based cuprates. However, the giant JT distortion stabilizes the CuC>4 plaquette with a 
big ex d x 2_ y 2 hole ground state which can be a proper starting orbital for an in-plane 
hole transfer b lg — > b lg in the CuC>2 planes of the quasi-2D cuprates. As a result we 
arrive at a disproportionation 2 CuO^ - — > CuOj~+Cu04~ with formation of a spin and 
orbital singlet electron center CuO^ - (analog of the Cu + ion) and a hole center CuO^ - 
(analog of the Cu 3+ ion) in a Zhang-Rice 1 Ai g ground state. Full disproportionation 
within a CuC>2 plane results in a system of composite two-electron spin-singlet (S-type) 
bosons moving on the spin-0 lattice formed by hole CuO^ - centers. We arrive at a 
EHBL system which is the most favorable one for the high-T c superconductivity. A 
comprehensive argumentation of the disproportionation scenario for quasi-2D cuprates 
is reported in Ref. [5j. 

At first sight the situation similar to cuprates can be realized in a hypotetical 
quasi-2D nickellate with a system of well isolated Ni 3+ C>2 planes provided a low-spin t\ g e g 
configuration of the NiO|~ center with the b lg oc d x 2_ y 2 electron ground state and a spin- 
singlet ZR-type configuration of the e 2 -subshell of the electron NiO^ - center. However, 
the latter condition seems unlikely to be realized in practice therefore the perspectives 
of the cuprate-like superconductivity in nickellates look like very illusive |57j . 

In the above we addressed "optimal" systems with bare octahedral <i 4,7 ' 9 complexes. 
What about predicted "optimal" d 1 and d 6 systems with tetrahedral local symmetry? 
At present we have no literature data about such d 1 systems whereas the systems with 
tetrahedral d 6 , or more precisely, high-spin Fe 2+ complexes number tens of the so-called 
iron pnictides ("FePn", where Pn is As or P) and iron chalcogenides ("FeCh", where Ch 
includes S, Se, and Te) |58j. Actually all these compounds are either superconductors or 
parent systems which turn into superconducting phase with T c values up to 56 K under 
nonisovalent doping and external or "chemical" pressure. 

All the families of iron- containing superconductors have two-dimensional planes 
of FePn4 or FeCli4 tetrahedra. The quasi-2D structure makes the iron compounds in 
common with superconducting cuprates, however, the FePn/Ch superconductors are 
fundamentally different from the cuprates in several points. First of all it concerns the 
coexistence of superconductivity with magnetism demonstrated at an unprecedented 
atomic scale. 

Unlike the case of the high-T c cuprates, electrons in FePn/Ch layers are always 
itinerant and there is no Mott insulating state in the electronic phase diagram though 
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most of the FePn/Ch systems start with an AFM order as the parent phase. The 
magnetic order in FePn/Ch layers is unconventional one, it is most likely a spin-density 
wave with strongly reduced mean values of a local effective spin. 

All the main features of the superconducting and normal state for the FePn/Ch 
compounds agree with our predictions of the disproportionation within the FePn/Ch 
layers accompanied by formation of EHBL system with a topological phase separation, 
e.g. formation of a multi-bubble system. However, in contrast with the spinless scenario 
for cuprates in the FePn/ Ch layers both composite electron e 2 g bosons and the hole lattice 
have nonzero spins, s = l and S = 5/2, respectively. It means that the charge ordered 
bubbles jointly are magnetic bubbles with a complex nanoscopically inhomogeneous 
spin structure. Accordingly, the FePn/Ch layers can be considered as a system of 
the bubble carriers which carry both charge, "local superconductivity", and a spin 
density. Moreover, the bubble formation is accompanied by strong nanoscale structural 
distortions which can be detected by pair distribution function analysis |59j . 




In any case, finding the high-T c superconductivity in the FePn/Ch compounds 
with the tetrahedral coordination of the iron Fe 2+ (3<i 6 ) ions in HS state and its 
coexistence with an unconventional magnetism can be a key argument supporting the 
disproportionation scenario |60] . Puzzlingly, our model implies the superconducting 
carriers in FePn/Ch compounds are composed of the e g rather than t 2 g electrons as 
predicted by one-electron band models. Both cuprates and ferro-pnictides/chalcogenides 
are believed to be "bad" EHBL metals. 

6. Conclusion 

We analyzed d-d disproportionation phenomena in materials with different 3d ions. 
We have shown that optimal conditions for the disproportionation driven high-T c 
superconductivity are realized only for several d n configurations. These are HS-ci 4 , 
LS-d 7 , d 9 configurations given octahedral crystal field, and d 1 , HS-<i 6 configurations 
given tetrahedral crystal field. Interestingly, all these bare systems are characterized 
by E-type orbital degeneracy, i.e. these are prone to a strong Jahn- Teller effect. In all 
the instances the disproportionation reaction lifts the bare orbital degeneracy, that is it 
has a peculiar " anti- Jahn- Teller" character. All the optimal systems are characterized 
by a strong suppression of the one-particle transport, while after disproportionation 
these are characterized by an effective two-particle transport. Anyhow the optimal 
disproportionation schemes (see Fig. |3] and Table I) imply a participation of a spin degree 
of freedom. Proposed superconductivity in the d A and d 6 systems at variance with d 1 , d 7 , 
and d 9 systems is characterized by unavoidable coexistence of the spin-triplet composite 
bosons and a magnetic lattice. Only quasi-2D cuprates in a "negative-?/" regime seem 
to represent a system of spin-singlet composite bosons moving on the spinless lattice. 
Their magnetic response could be explained by purely oxygen orbital magnetism if we 
take into account deviations from the simple Zhang- Rice model [291 US]- Interestingly, 
our "spinless" scenario for cuprates clearly opposes to a spin/magnetic mechanism to 
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be currently a leading contender for the superconducting mechanism in the high T c 
cuprates. 

We argue that unconventional superconductivity observed in iron-based layered 
pnictides and chalcogenides with tetrahedrally coordinated Fe 2+ ions can be a key 
argument supporting the disproportionation scenario is at work in these compounds. It 
is worth noting that the potential for high-T c appears to be the greatest for undistorted 
FePn 4 tetrahedra |61j. Furthermore, observation of the high-T c superconductivity in 
these quasi-2D systems with undesirable spin degrees of freedom both for the composite 
boson and the lattice gives some hope to finding of superconductivity in artificial 2D 
superlattices of all the "optimal" configurations. 

We did not concern here many important points of the disproportionation 
scenario and some comparisons with other model approaches, particularly with the 
bipolaronic theory by Alexandrov |41] . Nevertheless the model approach suggested is 
believed to provide a conceptual framework for an in-depth understanding on equal 
footing of physics of very different strongly correlated 3d systems such as cuprates, 
manganites, nickellates, ferro-pnictides/chalcogenides, and other systems with a charge 
transfer instability and/or mixed valence. In particular, our paper provides a clear 
answer to the question, "what is so special about Cu in cuprates or Fe in the 
Fe-based superconductors?" and why, e.g., BaMn 2 As2 which with respect to many 
properties is situated between BaFe2As2 and the high-T c cuprates is not a high-T c 
superconductor [61] . 
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